The critical perturbation parameter estimate 
for the transition from regularity to chaos 
in quantum systems 

V. E. Bunakov, I. B. Ivanov, R. B. Panin * 
Petersburg Nuclear Physics Institute, 188300 Gatchina, Russia 



Abstract 

In this paper we continue to develop our approach to the chaoticity properties of the quantum 
Hamiltonian systems. Our earher suggested chaoticity criterion characterizes the initial symmetry 
breaking and the destruction of the corresponding integrals of motion in a perturbed system, which 
causes the system's chaotisation. Transition from regularity to chaos in quantum systems occurs at 
a certain critical value of the perturbation parameter. In our previous papers we had to diagonalize 
the perturbed system's Hamiltonian matrix in order to estimate this parameter. In the present 
paper we demonstrate that the critical perturbation parameter for the transition from regularity 
to chaos can be estimated in the framework of the first order perturbation theory. The values of 
thus obtained critical parameter are in good agreement with the results of our previous precise 
calculations for Hennon-Heiles Hamiltonian and diamagnetic Kepler problem. 

1 Introduction 

For last decades the investigations of quantum chaos have been extensively carried out but 
this field remains under hot discussions. One part of researchers believe that quantum 
chaos doesn't exist and at best should be studied in the semiclassical approximation. It is 
beheved that the quantum analogue of the classical system should have properties refiecting 
regularity or chaoticity of classical trajectories. The law of level spacing distribution is 
considered to be one of such properties. It is also believed that the quantum analogue of 
the classically chaotic system obeys Wigner level spacing law, while Poissonian law holds for 
regular systems. However many authors (including us [1]), pointed out the incompleteness 
and crudeness of this criterion of quantum chaoticity. 

In this paper we continue to develop our approach [2-6] to the chaotic properties of the 
quantum Hamiltonian systems. Our main point is connection between the symmetry prop- 
erties of a system and its regularity or chaoticity. We had demonstrated that our earlier 
suggested chaoticity parameter characterizes the initial symmetry breaking and destruction 
of the corresponding integrals of motion in a perturbed system, which leads to chaotisation. 
As it was shown in [3], in the semiclassical limit the chaoticity parameter transforms into 
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Lyapunov's exponent of the corresponding classical motion. We had also shown that ac- 
cording to our chaoticity criterion the transition from regularity to chaos in classical and 
quantum systems occurs at the same critical perturbation parameter. 

In our previous papers we had to diagonalize the perturbed system's Hamiltonian matrix 
in order to estimate this perturbation parameter. In this paper we shall show that it is 
possible to estimate it in the framework of the perturbation theory. In part 2 we reconsider 
the definition of our chaoticity parameter which is used for the investigation of the transition 
from regularity to chaos . Part 3 is devoted to the analytical estimate of the critical parameter 
in the framework of the first order perturbation theory. The comparison of our estimate with 
the precise calculations for Hennon-Heiles Hamiltonian and diamagnetic Kepler problem is 
done in part 4. 

2 Destruction of quantum numbers 
and the chaoticity criterion 

Consider a stationary quantum system with Hamiltonian as a sum of Hamiltonian Hq of 
an integrable system and perturbation W , which destroys the symmetries of Hq in such a 
way (see e.g. [3-4]) that the number M of the independent first integrals of motion becomes 
less than the number of the system's degrees of freedom: 

H = Ho + \V. (1) 

Eigenfunctions {(pa} of the unperturbed integrable Hamiltonian Hq are common for some 
complete set of independent mutually commuting operators {Jp}p^i {N - the number of 
degrees of freedom of the Hamiltonian Hq). Eigenstates ipi of the full Hamiltonian may be 
expanded in eigenfunctions 0^ of the unperturbed Hamiltonian Hq. 

V'.=EC^^<^a = Ec?0a- (2) 
a a 

Let us consider the probability Pa{Ei) to find the basis state (p^ in the state i/jt with 
energy Ei, which is equal to the squared absolute value of the corresponding coefficient in 
(H) and define the energy width V^^^ of Pa{Ei) distribution: the minimal energy interval 
for which the sum of probabilities Pa{Ei) is larger or equal to 0.5. Thus defined F^^^ is the 
energy spreading width of basis state (pa- The spectrum of Hq may be degenerate and then 
the irreducible representations of the symmetry group of Hq consist of several basis functions 
which belong to one energy level (shell). We want to find a parameter characterizing the 
measure of initial symmetry breaking of Hq under the infiuence of the perturbation V. It's 
clear that such symmetry breaking results only due to significant mixing between functions 
from different irreducible representations. The mixing of states within one shell doesn't 
change their symmetry. While the spreading width V^^^ is smaller then the distance Dq be- 
tween the neighboring levels in the spectrum of Hq we can distinguish "localization domain" 
(in energy) of one set of basis states from "localization domain" of another one. When the 
spreading width exceeds Dq we start loosing the "traces" of basis functions in the spectrum 
of H and can't even approximately compare states ipi with irreducible representation of 
symmetry group Hq. Thus a parameter 

as" = T^/Dq (3) 
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is the natural measure of symmetry breaking. When the parameter se" exceeds unity the 
symmetry of the Hamiltonian Hq disappears. Such a value of perturbation is accompanied by 
disappearance of the initial selection rules, the levels are distributed approximately uniformly 
(level repulsion) and the level spacing distribution approaches Wigner's law. One can say 
that the transition from regularity to chaoticity has taken place in quantum system and ae" 
may be considered as the parameter of chaoticity. 

The spreading width r^^^ depends on the number a of the basis state of Hamiltonian Hq, 
i.e. on its quantum numbers. The classical analogy to this is the dependence of the invariant 
torus stability on the corresponding values of integrals of motion. It is clear that in order 
to obtain the global chaoticity characteristic in quantum case it is necessary to average F^^^ 
over the basis states (pa belonging to the same irreducible representation (shell). Thus we 
should determine the energy width F^pr- of the following distribution 



where T/v is a degenerate level subspace with main quantum number N. The averaged 
parameter F^pr unlike the local F^^^ has an important feature of invariance with respect to 
the choice of the basis for the integrable Hamiltonian Hq [1]. From the theoretical point of 
view the above chaoticity parameter has one more useful property. As it was shown in [3], 
in the semiclassical limit T^^^/h transforms into Lyapunov's exponent of the corresponding 
classical motion. Thus we see that parameter ?e{X,E) = Tspr/DQ "measures" destruction 
(fragmentation) of the irreducible representations of the basis and, hence of the Casimir 
operator (the main quantum number), which is the approximate integral of motion for small 
perturbations. 

3 Analytical estimate of the critical 
perturbation parameter 

The analytical estimate can be performed in the following way: Firstly, the expansion coef- 
ficients d^) at small parameter A can be computed in the framework of usual perturbation 
theory. Secondly, since it is inconvenient to work with the distribution width P^^Ei) ana- 
lytically, we shall use the equivalence of our chaoticity criterion based on and the criterion 
of Hose- Taylor approach. 

In a series of publications Hose and Taylor (see, e.g. [7,8]) in the framework of the the- 
ory of effective Hamiltonians suggested the criterion of their existence together with their 
integrals of motion. According to this criterion, one can construct a convergent sequence 
of approximations to the effective hamiltonian provided that a square projections of a per- 
turbed wave function on model space is larger than 0.5. Operators commuting with effective 
hamiltonian will be approximate integrals of motions of the problem. Thus, if a square of a 
projection of some states ipi of the Hamiltonian H on space the T/v of degenerate the level of 
Hq are larger than 0.5, then the principal quantum number must be an approximate integral 
of motion for these states (in this range of energies). Therefore the average square of a 
projection of perturbed wave functions on subspace T/v with fixed value of main quantum 
number should reach value 0.5 approximately at same critical perturbation when 33 = 1 [5]. 

The principal problem in using the perturbation theory for estimating of squares of pro- 
jections is the degeneracy of basis states and, hence, the necessity to solve a secular equation. 




(4) 
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which prevents to obtain the analytical expressions. In order to bypass this complexity we 
shall act as follows. We shall construct such a quantity, which on the one hand will help us 
to estimate the critical perturbation parameter, and on the other will be is invariant with 
respect to block unitary transformations inside subspaces T^- This property will allow us 
in evaluating this quantity to utilize not only the states obtained from the secular equation, 
but also any others, generated from them by the block unitary transformations, (including 
the initial basis states 0^). In this way we shall not neeed to solve the secular equation. 

Let us estimate the average square of projection W{\, N) of the perturbed wave functions 
not on a subspace 7/v, but on its orthogonal adjoint HgT^. Further, let us assume, that we 
already have solved the secular equation in all subspaces TJv and found correct functions (p^^ 
in the zero approximation. Then in the first order approximation the expansion coefficients 
of perturbed states i on basis a belonging to other shells will be: 

(5) 



where Vai = (f)''^^*V(j)f''^ and E^^ is the energy of states in a zero approximation. The square 
of projection of a perturbed state ipi on orthogonal adjoint to T^r will be 



IK,: 
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Now we average the previous expression over all dimTiq states ipi (for small perturbations 
every state tpi can be approximately referred to some 7/v): 

With increasing perturbation the parameter W{\^ N) will gradually grow (the basis states 
are fragmented over other shells) and at some critical value will reach 0.5. Thus a required 
critical perturbation parameter can be found by solving the equation 

M/(A,iV) = 0.5 (7) 

with respect to A or A^. 

Let us prove now the invariance of W{\, N) with respect to arbitrary block unitary 
transformations of the basis, which mixes the functions 0^°-' only inside the irreducible rep- 
resentations 7/v ( the functions from different T/v are not mixed). It is possible to present 
the sum of the states a ^ T/vas 

E = E E • (8) 

o^Tjv n^^N a&r„ 



Therefore (H) can be written as 



W{\N) = -^Y.Wn, (9) 



IK 
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= T. H (10) 
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Now consider what happens with the value Wn under the block unitary transformation Ui, 
of the zero approximation basis wave functions 

#^0r = Ef^.4°^- (11) 

Since the energies E^^^ and E^^ are the same for all states i & Tn and a & %, and do not 
vary under the transformation the denominator in (|lO|) we can be taken out of the sum. 
Utilizing the symmetry in labels i and a we can prove the invariance of Wn with respect to 
their exchange. Indeed, 

E E 0r*v^Ef^.4°^Ef^:.(4°^*^€^) = E0f*^E(Ef^M^ = 
E€^*^E0r(<^r^<^i°^) = EE iv^«.r = (4°^ - E^^^m. 

a fi a fi 

Thus, making transformation Ub at first inside the subspaces 7/v, and then inside we 
obtain the invariance of Wn with respect to Lfb. Since each term Wn in the sum (P) is 
invariant, all the sum W{X, N) will be also invariant with respect to unitary transformations. 

This property allows us to bypass the necessity to solve the secular equation. We can use 
in the evaluation of W{X,N) in Eq.(P) the initial basis functions (pa and energy Ej^\ and 
find the critical perturbation parameter value from the equation (^. 



4 Two examples 

In this part we shall compare the analytical estimates of the critical perturbation parameter 
with the precise calculations for a nonlinear Hennon-Heiles Hamiltonian and diamagnetic 
Kepler problem. 

Let us consider a known Hennon-Heiles system with hamiltonian: 

H{q,p) = lipl + ql) + \{pI + ql) + \{qlq2 - ql/S). (12) 

Introducing operators of creation and annihilation al = -^{qk + ipk), dk = '^{qk — ipk), k = 
1, 2 we construct, as usual, the two-dimensial Cartesian oscillator basis 

0^ = = ^^(4)"n4)"¥o, (13) 

Eni,n2 = En = Hrii + ^2 + 1). 



It is easy to write out matrix elements of the Hamiltonian ( [121) in the basis (0) and to find 
through a diagonalization the system's eigenergies and eigenstates. In the calculations of 
bound states (with energies E < 1/(6A^)) the value A = 1 and h = 0.01 was fixed and the 



basis ([T3|) was used, which included 496 of states (30 shells). We investigated average (over 
states with approximately equal energy) square of a projection of perturbed wave functions 
on subspace TCqT^. First we carried out precise calculations (see points on fig. |I]), and then 
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Figure 1: Perturbation theory (line) and precise calculation (points) of average square of a 
projection of wave functions on a subspace TiQTjq dependence as functions of energy in the 
Hennon-Heiles system. 

determined the value W{\, N) from Eq.(^) (line on Fig. |l|). According to precise calculation 
we achieve critical value 0.5 at an energy E ^ 0.105, that practically coincides with earlier 
obtained [4] value with the criterion ffi = l(£'^0.11). As one should to expect, the 
perturbation theory works well at small perturbation parameter (up to energies E k, 0.08) 
and as a result yields the critical value of energy E k, 0.084. 

The same test of the critical perturbation parameter estimate was carried out for dia- 
magnetic Kepler problem with Hamiltonian 

H = p^/2m - (?lr + ujI^ + \^mJ'{x^ + y^). (14) 

Here frequency lo = eB/2mc, B — magnetic field, directional along an axes z. Instead of the 
usual energy E we used the standard dimensionless energy e = where = huj/n (7^ 

— Rydberg constant). We diagonalized a matrix of the Hamiltonian in the Coulomb parabol- 
ical basis with quantum numbers {ni,n2,m) and energies of basis states E^ = —l/2n'^, 
n = ni + n2 + \m\ + 1. For the investigation of a fragmentation of basis functions with m = 
from a shell with the principal quantum number 10 we used the basis of the first 20 shells. 
We studied the dependence of the average square of a projection of the perturbed wave 
functions on a subspace 7i 7/v as a function of the scaled energy. As in the Hennon-Heiles 
case, we fixed a subspace Tiq with a main quantum number = 10 and varied the values 
7. As an averaging range we selected the part of the spectrum, whose states {dimTio = 10) 
would have the maximum square of a projection on Tiq. Fig. ^ shows the analytic result 
W{X, N) (line) and the numerically clculted (points) dependence of the precise wave func- 
tions' average square projections on a subspace H Tiq as functions of the scaled energy e. 
As well as in the previous case we have obtained the quite good agreement of the theoretical 
critical parameter estimate [ecr = —0.54) with precise calculation (ecr = —0.47), which is 
also in correspondence with earlier results [6] obtained using = 1 criterion {ecr = —0.45). 
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Figure 2: Perturbation theory (line) and precise calculation (points) of average square of a 
projection of wave functions on a subspace HqTiq dependence as functions of scaled energy 
for a diamagnetic Kepler problem. 

5 Conclusion 

In this paper we continue to develop our approach [2-4] to the chaotic properties of the 
quantum Hamiltonian systems. Our main point is the connection between the symmetry 
properties of a system and its regularity or chaoticity. We show that the earlier suggested 
chaoticity parameter ^(A, E) characterizes the initial symmetry breaking and destruction of 
the corresponding integrals of motion in a perturbed system, which leads to chaotisation. The 
value of ce is equal to the ratio of the average spreading width of the basis wavefunctions over 
the perturbed ones to the average unperturbed level spacings. In order to estimate the critical 
perturbation parameter for the transition from regularity to chaos we had to diagonalize the 
perturbed system's Hamiltonian matrix. In the present paper we have shown that the critical 
perturbation parameter may be estimated in the framework of the first order perturbation 
theory. We construct a quantity which is invariant with respect to unitary transformations 
of subspaces T/v of basis degenerated levels Hq and indicate the measure of system chaoticity. 
Thus there is no need to deal with secular equation and the approximate value of the critical 
perturbation parameter might be obtained from the relation which contains matrix elements 
of perturbation V in the initial basis and energies of basis states. The values of thus obtained 
critical parameter are in good agreement with the results of our previous precise calculations 
for Hennon-Heiles Hamiltonian and diamagnetic Kepler problem. 

One of the authors (IBI) is indebted to Prof. Zikiki and to the Organizing Committee of 
V. Gribov's Foundation for their support. We also thank RFFI (grant N 00-15-96610). 
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